function [alpha_abs,alpha_abs_Int]= AbsCross(detuning,inten)

% [alpha_abs,alpha_abs_Int]= AbsCross(detuning,inten)
% Calculates the absorption cross section for the Rb D2 transition
% 
% Outputs:
%       alpha_abs = absorption cross section 
%       alpha_abs_Int = absorption cross section, taking into account possible intensity saturation

cons = Constants();

alpha_abs = (3*cons.lambdaD2^2/(2*pi))/(1+(2*detuning/cons.Gamma)^2); % abs. cross section
alpha_abs_Int = (3*cons.lambdaD2^2/(2*pi))/(1+(inten/cons.Isat_lin)+(2*detuning/cons.Gamma)^2); % abs. cross section with intensity saturation taken into account